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In this paper we are concerned with the simplest normal form computation of the systems 

x = 2xf(x, y 2 + z 2 ), y = z + yf(x, y 2 + z 2 ), z = -y + zf(x, y 2 + z 2 ), (0.1) 

where / is a formal function with real coefficients and without any constant term. These 
are the classical normal forms of a larger family of systems with Hopf-Zero singularity. 
Indeed, these are defined such that this family would be a Lie subalgebra for the space of all 
classical normal form vector fields with Hopf-Zero singularity. The simplest normal forms 
and simplest orbital normal forms of this family with non-zero quadratic part are computed. 
We also obtain the simplest parametric normal form of any non- degenerate perturbation of 
this family within the Lie subalgebra. The symmetry group of the simplest normal forms are 
also discussed. This is a part of our results in decomposing the normal forms of Hopf-Zero 
singular systems into systems with a first integral and nonconservative systems. 

Keywords: Parametric normal form; Hopf-Zero singularity; Eulerian vector fields. 



1 Introduction 



Normal form theory is one of the most effective tools for the local study of nonlinear dynam- 
ical systems. The basic idea is to use permissible transformations and obtain a simplified 
vector field. Transformations are permissible that preserve certain dynamical features of 
the original system. The space of all permissible transformations form a group and acts 
on vector fields like an action of a group on a vector space. Consider a set of vector fields 
generated by the group acting on a given vector field. Then, the infinite level normal form 
of the vector field is to find a unique representative from this set. Thereby, the computation 
of infinite level normal forms is an important tool for classification of vector fields. The 
uniqueness of a normal form computation is determined by the (infinite) level of normal 
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form through an specifically chosen normal form style and costyle. The level of a normal 
form assesses the remaining spectral data at our disposal in the permissible transformation 
space for further simplification of the system while a normal form style and costyle makes 
a unique choice for the normal form vector field in each level of normal form computation, 
see e.g., [HI EH [25]. Therefore, the infinite level normal form is sometimes called the sim- 
plest normal form or unique normal form when a normal form style have already been fixed. 
When a system has some symmetries, it is important that its normal form would preserve the 
symmetries. Although there are research results on the simplest normal forms of symmetric 
systems but they are considerably less than the existing results on normal forms without 
symmetry. There are several reasons for this. The first difficulty is to recognize the symme- 
tries and then, to find the group of transformations preserving the symmetries. Therefore, 
one is also concerned with the space of the symmetric vector fields invariant under the group 
action. This needs a good understanding on the algebraic interactions of the symmetric and 
nonsymmetric vector fields with the transformation groups. Once all these are successfully 
accomplished, in most cases the normal form computation is more difficult in systems with 
symmetry than in systems without symmetry. 

Normal form decomposition of a nonsymmetric vector field into two symmetric vector 
fields can have many important potential applications. For example Eulerian and Hamilto- 
nian vector fields are two important families of vector fields. Therefore, the study of their 
dynamics is important. Normal form decomposition of arbitrary vector fields into Eulerian 
(nonconservative or dissipative) and Hamiltonian (conservative) vector fields are also impor- 
tant in both theory and applications. Wiggins [30| Chapter 33] remarks that transforming 
a system into an integrable Hamiltonian system plus a nonconservative perturbation facil- 
itates "a wealth of techniques for the global analysis of nonlinear dynamical systems such 
as Melnikov theory, perturbation theory for normally hyperbolic invariant manifolds, and 
Kolmogorov, Arnold, and Moser (KAM) theory" . Furthermore, Palacian [22] indicated that 
such kind of decomposition can be used in developing an ODE solver. For example using the 
obtained scalar function as the integral of a piece of the vector field in an ODE solver, can 
enhance its efficiency These signify the importance of developing methods for such kind of 
decomposition. Recently, some researchers have paid attention to this theory and have made 
important contributions to the subject, see e.g., [3, 4, 7, 8, 16, 21-25, 27]. Furthermore, 
the dynamics of a system may sometimes be analyzed by dynamics of its components. For 
example, an advection-diffusion model consists of advection terms and diffusion terms and 
its dynamics is well understood as a combination/competition of the dynamics associated 
with the advection and diffusion terms, see e.g., [HI Section 5.5]. Therefore, the study of 
each components of a decomposed vector field may help to better understand the dynamics 
of the full system as a combined dynamics or as a competing behavior between its compo- 
nents. Thus, it is important to individually deal with the cases that the vector field is a 
quasi- Eulerian, or is a conservative vector field, see [10]. In this paper a vector field is called 
conservative when it has a first integral and is called nonconservative when it does not have 
any first integral. 

This paper deals with the nonconservative family of our upcoming results on such de- 
composition for Hopf-Zero singular vector fields. Systems with Hopf-Zero singularity are 
important in applications. There are several important research results on the simplest nor- 
mal forms of Hopf-Zero singularity, see [H El El |27]. However, there does not seem to exist 
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any result on such kind of normal form decomposition for Hopf-Zero singularity. Although 
the Hamiltonian vector fields (symplectic structures) require an even dimensionality, the 
decomposition idea still can work. Indeed we instead work with conservative and noncon- 
servative family of vector fields, see [TT]. In this paper we are concerned with the study of 
nonconservative family. Indeed, we consider the family of Hopf Zero singularities given by 

d d 

V = %- y d-z +V ^ (L1) 

where 

and / is a formal power series function and without constant terms. Any result in this 
direction also helps in a better understanding of the more complicated algebraic structures 
involved in the three dimensional nilpotent singularities. In this paper the coefficients are 
taken from a field F of characteristic zero, see [25] where the coefficients are taken from a 
local ring containing the rational numbers. Since E is a quasi-Eulerian vector field, is 
called a formal quasi-Eulerian vector field, that is a nonconservative vector field. The reader 
should note that this family of vector fields are indeed the classical (first level) normal form 
of a very much larger family of vector fields. This, however, is independent from [10J and 
basically refers to an essentially separate problem. Our normal form computation is to use 
a transformation group that preserves the Eulerian structure and it would only simplify 
the function /. However, our results in [10] work with the bigger group, that is, the set of 
all near-identity transformations and may create conservative terms into the normal form 
system. Indeed, the quasi-Eulerian vector field E is defined such that its generated family 
would form a Lie algebra. 

This paper is organized as follows. Section [2] introduces a family of nonconservative 
vector fields that is a Lie algebra and presents a general notations for normal form compu- 
tations. The simplest normal forms and orbital normal forms of this family of vector fields 
with non-zero quadratic part are obtained in Sections |3] and EJ respectively. In Section [5] we 
prove that under certain technical conditions multi-parameter nonconservative perturbations 
of the vector field v given by Equation (11.11) can be transformed into 

^ (oo) = zj- - y-^ + g(x, y 2 + z\ fiE, (1.3) 

where g(x, R, /x) = x + ffiR + £[=o + Pr Rr for a r > 2. 

2 Nonconservative vector fields and notations 

This section presents the necessary Lie algebra structure and notations that we need for 
normal form computations in the following sections. 

Definition 2.1. Define 

e[ : = v + z 2 r 1 ^ + \*\y 2 + - 2 )%^ + \Av 2 + (o < i < m-i) 
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Since E{ = x l (y 2 + z 2 ) k ~ l E and E is a quasi- Eulerian vector field, we call E k a quasi- Eulerian 
vector fields. Then define 

S := spanj I Pi ^ F > 1 + k > h k > l}. (2.2) 

Let 9q := z-^ — y-j^. Therefore, the vector field given by Equation (11.11) can be represented 
by an element of 0° + $ . 

The following theorem follows from a straightforward computation. 

Theorem 2.2. The vector space {clOq + v\a G F, v G $} is a Lie algebra with the following 
structure constants: 



d+m 
k+n i 



[ElK] = (n-k)E, 

[E l k ,e° ] = 0, 

for any I, k,m,n G N . 

Proof. The space $ is a Lie algebra by the Lie bracket [v, w] := vw—wv for any v,w G <S . The 
structure constants follows from the fact x l (y 2 + z 2 ) k ~ l is an i^-eigenfunction with eigenvalue 
k. □ 

Note that one may consider the ring generated by polynomials x l (y 2 + z 2 ) k ~ l . Then, it 
forms a Lie algebra; the Lie bracket of two eigenf unctions is given by their product mul- 
tiplied with the substraction of their eigenvalues. Thus, normal form computations can 
be implemented on a computer algebra system using scalar valued functions rather than 
using vector fields. This enhances the efficiency of the computer program. The space 
{aOg + v\a G F, v G #} is indeed a Lie subalgebra of all the classical normal form of 
Hopf-Zero singularity, see e.g., [T| |5| |6| ITT]. 

Theorem 2.3. There does not exist any first integral for any v G <f , i.e., any v G $ is a 
non- conservative vector field. 

Proof. Suppose that / = J2 m k a rn>n ^x m y n z k is a first integral for gE G S where the 
function g is nonzero. Then, f x x + f y y + f z z = implies a m , n ,fc = for all m,n,k G No, i.e., 
f — 0. This completes the proof. □ 

Now we briefly present a general framework on how to obtain the simplest (orbital and 
parametric) normal form of a vector field, for more detailed information, terminology and 
background material see [14]. Let $ = ^ ^ e a g ra ded Lie algebra, i.e., [vi, Vj] G for 
any v\ G Si and Vj G Sj. The sub- indices i for spaces (not element indices, of course) in this 
paper refer to the homogenous spaces of grade i according to the associated grading. Recall 
that a style of a normal form is a rule to choose a unique complement space to a subspace 
of $ . Assume a normal form style is given and A = A i s the graded vector space of all 
generators associated with the permissible transformations (in [H] this is denoted by A ). 
The space A acts linearly on S . This action is denoted by * and it preserves the grading 
structures, i.e., At* Sj C £ i+ j for any i,j. 
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Remark 2.4. The action * is the initially linear map of the associated near-identity per- 
missible transformations generated by elements of A, see [H] for more detailed discussion. 

The space A and its associated action * are different for the normal form, orbital normal 
form and parametric normal form of vector fields. Therefore, we will explicitly define these 
accordingly in the following three sections. Let v = ^°^ Vi and inductively define 

d nA . An by d n,i( Yn ) := Y n * v for Y n G A n , (2.3) 

and 

d n ' k : A n x ker{d n ~ 1 ' k - 1 ) ->■ S n (for any jfc < n), (2.4) 
given by d n ' k (Y n , Y n _ u Y n _ k+1 ) := Y^Zl Yn-% * «», where 

(y n _ x , y_ 2 , . . . , y n _ fc+1 ) g ker(d"- 1 - fc - 1 ). 

For any jfc > n, let d n ' fe := <f« n . Note that (Y n _ u Y n _ 2 , . . . , Y n _ k+1 ) G ker(d n " 1 ' /c - 1 ) if and only 
if 

3 

^2 Yn+i-j-i * Vj-i = (2.5) 

i=Q 

for any j = 0, 1, . . . , k, see also [29] . Let TZ n,k := Im (d n,k ). Then, there exist the complement 
subspaces C n ' k such that 

K n ' k © C n ' k = <§ n (2.6) 

follows the normal form style. We call d n,k the &-th level state map. If Vi G C n,k for any 
natural number n, we say that f = f j is in the fc-th level extended partial normal form (see 
[T8] for the original definition and ideas) and when Vi G C l,i for any natural number i, v is 
said to be in the infinite level (orbital or parametric) normal form. Notice that the domain of 
the map d n ' k has ^-components. However, the first component does not play any role in our 
normal form computation, since eg is the zero grade of the system and we have [$, 0[j] = 
in Section [HI and in Sections [4] and [5] the linear part 0q is first removed before normal form 
computation. Thus, we omit the first component and we only have (k — l)-components in 
the domain of the map d n,k , i.e., d n ' k (Y n _ 1 , . . . , Y n _ k+1 ) = Yli=i Y n -i * Vi- 

Theorem 2.5. There always exist a sequence of permissible transformations transforming a 
vector fields into its k-th (or infinite) level extended partial normal form. 

Proof. The proof follows from [14, Theorem 4.4]. □ 

In this paper, we also discuss the symmetry group of the infinite level normal forms. The 
symmetry group of a vector field is a subgroup of the permissible transformations such that 
the infinite level normal form invariant under the symmetry group acting on v^°°' . This 

symmetry group is generated by the infinite level co- normal form space. The infinite level co- 
normal form space is a subspace, say B, of the space of permissible transformation generators 
A in which B * = and can be computed through the normal form computation. 
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3 The simplest normal forms 

In this section we obtain the simplest normal form of vector fields given by Equation (II. ip 
with nonzero quadratic terms, i.e., 

V W ■- 0° + p\E{ + ■ ■ ■ E 0° + S, (3.1) 

where j3\ ^ 0. Define the grading function 5 for the space $ by 

8{E l k ) := 2k -I and 5(Q° ) = 0. 

Then, S(E l k ) > k. For simplicity of notation 5 or 5 a is used for all gradings in this paper. 

Lemma 3.1. There exists a formal near identity transformation that sends the vector field 
v given by Equation 113. 1\) into the second level extended partial normal form 

oo 

= 0° + P{E\ + 0°E$ + PlEl + &E\ + PlEl (3.2) 

k=2 

where (3[ is not changed. 

Proof. Define A := Y^i^i- Since [E l k ,E\\ = — (k — l)E l £\, a complement space to the 
image of d n ' 2 is span{i?2fc} for any n = 2k > 3, spanji?^} for n = 3, and span{i?f , E®} for 
n = 2. □ 

The reader should note that in this paper the notation (3 l k for coefficients are not changed 
as they are updated in the normal form computation. 

Theorem 3.2. There exists a natural number r and a formal near identity transformation 
sending v^ 2 ' given by equation h3.2\) into one of the following (2r + 2)-th level extended partial 
normal forms: 

oo 

v (2r +2) = Q o + plE{ + + m + fc E i + p o E o + ST fiE* (for r > 2) , (3.3) 

k=r+2 

or 

oo 

v (2r + 2) = Q o + p i E i + pa E o + m + m + m _ m)E o + + J- $El (3.4) 

k=r+2 

where f3i((3\ — /3i(3 2 ) ^ and r > 2. Furthermore, t/ 2r+2 ) = v^°°\ and for /3° ^ there is no 
non-trivial symmetry in the symmetry group for v^. 

Proof. Consider the vector field obtained in Lemma [3. 1[ Denote 

£o,k+i ■= (o, . . . , 0, - (3l(3 2 2 )El $El El o) 

and 

S 1M2 := (0, . . . , 0, - - $%))El (pi - ftPDEl ftEl E\) , 
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where the number of zeros is (k — 3)-times for any k > 3. Then, we have 

ker(<i 5 ' 5 ) = span{£ 0i5 }. 

In the case that $ - - ffiffi) ^ 0, we have 

d^{£ 0fi ) = (# - 0°M - m))El G W 6 < 6 and <?>\£ lfi ) = - - (3 2 )E 3 . 

Then, f3®£ Q + £i,6 £ kerrf 6,6 . Let r := 2, and hence the 6-th level extended partial normal 
form is given by Equation (13. 3p . 

For /3 2 ° = 0*(p\ - 131131), let r := min{£; | (3° k ^ for k > 2}. When /3° = for all k > 1, 
the choice r := 2 trivially completes the proof. Thus, we can assume that (3® ^ and r is 
finite. Therefore, ker d k,k = span{£ ,fc, £i,k} for any 6 < k < 2r+ 1, 

^+2 )2 , + 2 ( ^ 2r+2) = (r _ 1)(3 E + ^ and ^ + 2, 2 r +2(fi 2r+2) = _ (r _ l)^0 E +l . 

Thereby, £ r ° +1 G Im d 2r+2 ' 2r+2 and ker rf 2r+2 - 2r+2 = span{£ 2 r+ 2 }, where £ 2r+2 := /^£ , 2r+2 + 
£\ y 2r+2- Hence, ker<i fc,fc = span{£fc} for k > 2r + 2 is one dimensional. Identify with an 
element in S ', then converges to the infinite level normal form vector field in the filtration 
topology when k approaches infinity, see [13J. This proves that the simplest normal form 
computation does not produce any non-trivial symmetry for the normal form system. The 
proof is complete by Theorem 12.51 □ 

4 The simplest orbital normal form 

In this section, we obtain the simplest orbital normal form of v given by Equation (I3.ip . Let 
Z™ := x m y n ~ m for m < n,m,n G N, and Zq := 1. Then, 

M := span j a n Z nK G F, n ^ o} 

n>m>0 

is the space associated with time rescaling generators, i.e., ¥[[x,y}} = M © F. The time 
rescaling computation is governed by a module structure, see e.g., [T4] . Thus, the structure 
constants for the ring M and the ^-module <§ are given by 

ZnZ[:=Z^l and Z™E l k := Eft*. (4.1) 

Furthermore, let 

r := min{fc|/3£ ^ 0}. (4.2) 

A grading function is defined by S(Z l k ) := 2k — I and thus, M = Define Aj := Si x Mi 

for any i > 1, A := YjtLi A- Here A is a graded vector space according to [131 E] an d denotes 
the set of all formal sums $2i^i( u *> where G <^ and T; G The action of A on is 
given by (5, T) * v :— [S, v] + Tv for any (S, T) G A and v E S. 

Theorem 4.1. There exists a natural number r such that the vector field v G 0[j + S with 
nonzero quadratic part can be transformed into the simplest orbital normal form 
v (oo) = Q o + E i + + pO E o^ for any r > 2 

Furthermore, the symmetry group of the vector field v^°°' is not finitely generated. When 
f3® the symmetry group has no non-trivial symmetry from the group of permissible 
changes of state variables. 
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Proof. We first remove the linear part 0g from the system with a linear change of state 
variable [x, y, z] := exp{{Z-^ — Y-j^)t) [X, Y, Z], sending the new variables [x, y, z] := x-j^ + 
y-r^ + z^ into the old variables, see also [T7J the proof of Lemma 5.3.6]. Therefore, we may 
assume v := E\ + . . . G $ . Since 

d 2 ' 2 (0, Z{) = -$El G Im<2 2 ' 2 and <2 3 ' 2 (0, Z\) = - /3 2 2 /3 1 °) J B 2 1 G Imd 3 ' 2 , 

both terms of E 2 and i*^ are eliminated from the second level extended partial orbital 
normal form, i.e., — E\ + YlT=i ^k^k- Now let r : = min{k\(3 k ^ 0, k > 1}. (In case r 
is not finite, an arbitrary choice for r > 1 trivially completes the proof.) Denote &mk+2 : = 
(0, . . . , 0, Z^, ■j-j3^y-E'm) where the number of zeros is 2/c-times for k > 1. Then, 

<s = (0,0,21,^) e kerd 2m+2 > 3 and ^ +2 ^ +1 « r ) = ^/3 r °^ +r , (4.3) 
V (1 — m) J m — 1 

for any m > 1. By Theorem 12.51 there exist transformations that send the vector field t> into 
the infinite level orbital normal form u := E\ + /^.E^ + /3°£y G ^ for some r > 2, /3* G F. 
For any Z l k (I > 1, k > 2) there exists a nonzero state solution S £ y and a time solution 
T G J such that (T, S 1 ) * u = — Z k u. Therefore, the symmetry group for u is not finitely 
generated. Now the linear change of state variable [x,y,z] := expiiY-^ — Z-^)tj [X, Y, Z\ 
puts the linear part 0° back into the normal form system, i.e., Equation (14. 31) is 

obtained. The composition of this map with the symmetry group associated with u results 
in a symmetry group for v^°°\ The rest of the proof follows Theorem 13. 2[ □ 

5 Parametric normal form 

This section is concerned about parametric normal form of a non-degenerate multi-parameter 
nonconservative perturbation w (/i) (will be defined bellow) of the vector field v given by 
Equation (fLTft . see also PEJUS]. Define 

y := span{ ^/^ m i^ m |/t m G f}, (5.1) 

where the summation is over m G Ng, I + k > 0, k > I, := (//i, (/,%, . . . , fJ, q ), |m| = 
Sj=i m i — 0'/^oo = 0' anc ^ A 4 ™ := Z 1 ™ 1 ■■■lH q - We similarly define the parametric space 
^ such that F[[x, y, /i]] = © F, and the space ^ as the space of all formal power se- 
ries in terms of parameters \x G ¥ q without constant term, i.e., ¥[\p]] q = © ¥ q . The 
algebraic structures are naturally extended from nonparametric to the parametric cases, see 
e.g., [mm]. For m = (mi, . . . , m q ) G Nq and a = (oti, . . . , a 9 ) G N 9 , let am := £)j=i oijrrij. 
Then, define the grading function by 

5 a (E l k {i m ) :=2k-l + a-m, 5 a (Z l k [i m ) := 2k - I + a ■ m, and 5 a (// m ) := a • m. 

Define 



oo 

:= SA x x ^, A:=Y^ A, and (£, P,T)*v := + + [5, u] 

i=l 



M. Gazor, F. Mokhtari and J. A. Sanders 



Hopf-Zero singularity 9 



for any (S, P, T) G A. The perturbation of the system is considered to be within the space 
of i.e., w(n) G 9[j + y . We call any vector field 

w(fi)eQ° + y (5.2) 

a multi-parameter nonconservative perturbation of the vector field v given by Equation (11.1 p . 
when w(0) = v. Similar to the previous section, we remove the linear part 0$ from the vector 
field with a linear change of state variable. Since Wq = 0, the first component of the map 
d n ' k is omitted, e.g., the map 

<T' 2 : A n ^ -> S n , (5.3) 

is defined by d n ' 2 (Y n -\) = Y n -\ * V\ for any Y n G A n . Since f}\ ^ 0, through a time rescaling 
we may assume that (3\ Q = 1, i.e., w = Y^Li w i-> w i = anc ^ w i = S ^ m-^fc/^™ with the 
summation over /, fc, m G Nq, i = 2fc — i + a ■ m for any i > 1. 

Lemma 5.1. .Am/ multi-parameter nonconservative parametric perturbation of a vector field 
given by equation U.l\) can be transformed to the second level extended partial parametric 
normal form of given by 

w^ = El+ Yl 4° m ^V m (&°o = 0). (5.4) 

fc>0,|m|>0 

Proof. Let a :— (1, . . . , 1). Since 

ZVi = 0, [4, £7*] = (1 - k)E l + + \, Z°El = El Z\E\ = El and Z\E\ = E 2 2 , 
for any N > 2 we have 

lmd N ' 2 = span{£[// m | for any I > 1, Jfc > Z, where N = 2k — I + |m|}. (5.5) 

Thereby, ^T^' 2 = span{£^/x m | JV = 2fc + |m|} for any TV > 2. □ 

Consider the second level extended partial parametric normal form given by Equa- 
tion (15.41) . Let r := mm{k \ (3 k ^ 0, k > 1}. For the rest of this paper, we assume that there 
exists /?°o 7^ f° r & k > 1, i.e., r is finite. Define a (r + 1) x q matrix 

A := (fly), where := /3°_ 1)e . and ej denotes the j-th standard basis of ¥ q . (5.6) 

Then, we call a parametric vector field u>(/u) G Oq + a non-degenerate multi-parameter 
nonconservative perturbation of the vector field v if v = w(0), and rank(A) = r + 1 when 
w(/x) is transformed into given by Equation (15 ,4p . When this condition is satisfied, 
through an invertible linear change of parameters we can transform the vector field into 

r r 

w^ = E 1 1 +J2^ E° k + J2E^ k+1 + 4°m^V m + E ^V^V m , (5.7) 

fe=l fc=0 r->fc>0,|m|>2 fc>r,|m|>0 

i.e., aij — if i =fi j, and a,, = 1. Now we update the grading weight vector a by a; : = 
2(r — % + 1). Therefore, 5 a (.E^/ifc + i) = 2r for any < k < r. 
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Lemma 5.2. Let w(p) be a non- degenerate multi-parameter nonconservative perturbation 
of the vector field v governed by Equation U.l\) . Then, there exist a natural number r and a 
sequence of permissible transformations such that they transform w into its (2r + l)-th level 
extended partial parametric normal form 

r 

w Vr + i) = E i + p iE o + p rE o + J2E° k Vk + i (5.8) 

fc=0 

Furthermore, w^°°' = w^ 2r+1 \ The symmetry group of is not finitely generated, yet it 
has no non-trivial symmetry from the group of permissible changes of state variables. 

Proof. Without loss of generality we assume that w is transformed into given by Equa- 
tion E3 and (3% + 0. By Equation (TOD, we have E° m+r p m e Im d 2m+2r + a - m ' 2r+1 for any 
m > 1. Furthermore, the condition rank(A) = r + 1 implies that E®p m 6 Im c / 2fe + a - m > 2r + 1 for 
any |m| > 2,0 < k < r. Theorem 12.51 implies that there exist transformations such that 
can be transformed into y/ 2r+1 ) as desired. Since ^( 2r+1 ) does not have any term of grade 
higher than 2r, we have it/ 00 ) = w^ 2r+1 \ The uniqueness of together with an argument 
similar to the proof of Theorem 14.11 complete the proof. □ 

A linear change of state variables creates the linear part 0{j back into the infinite level 
normal form system. This together with Lemma 15.21 proves the following theorem. 

Theorem 5.3. Let the hypothesis of Lemma \5.S\ hold. Then, there is a sequences of permis- 
sible transformation transforming the system generated by w(fi) into the infinite level normal 
form system 

x = xg(x,p,n), p = ^pg(x,p,^), = 1, (5.9) 

where g(x, p, p) = x + (5ip 2 + J2l=o P 2 Vfc+i + PrP 2T ', A^O and r > 2. 

The following theorem sheds light to the dynamics of the infinite level parametric normal 
form system. 

Theorem 5.4. Consider the infinite level parametric normal form system in Theorem \5.3[ 
Then, the surface p 2 = c\x\, for any nonnegative constant c, is an invariant manifold of 
solutions for the system. Furthermore, the surface x = —fi\p 2 — J2k=o P^t^k+i ~ PrP 2T is an 
invariant manifold of limit cycles for the parametric system; except that it includes the origin 
as a fixed point. This implies that any point outside the later invariant surface, approaches 
a limit cycle or the origin through the surface p 2 = c\x\ for an appropriate constant c. 

Proof. By Equation (12. ip we have ^px = xp and thus, p 2 = c\x\ is an invariant set for the 
system. The rest of the proof follows from Theorem 15.31 □ 



Now we present an example and find its parametric normal form. 
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Example 5.5. Consider the system governed by 

x = 2x/ii + 2x 2 {a + fn+ /i 2 ) + 2x{y 2 + z 2 ){b + /i 2 ) + 2x 3 {c + fij + /i 2 2 ) + 2x(y 2 + z 2 ) 2 > 3 
+2ex 4 + 2x 2 ( ? / 2 + 2 2 )((i + ^) ) 

y = z + ynx +xy(a + ni + /i 2 ) + y(y 2 + z 2 )(b + /j, 2 ) + x 2 y(c + fi 2 2 ) + y(y 2 + z 2 ) 2 fi 3 

+ex 3 y + xy(y 2 + z 2 )(d + /i 2 ), 
i = -y + zfa + xz(a + /ii + /i 2 ) + ^(y 2 + z 2 )(b + /i 2 ) + x 2 ^(c + /i 2 2 ) + -2(?/ 2 + z 2 ) 2 ^ 

+ex 3 z + x^(y 2 + z 2 )(d + /i 2 ). 

where ab(da — cb) ^ 0. Then (using a Maple program) we obtain its parametric normal form 
given by 

x = 2xg(x,y 2 + z 2 ,/j,), 
V = z + yg(x,y 2 + z 2 ,fi), 
z = -y + zg(x,y 2 + z 2 ,/j,), 

where g(x, R, fi) := ^ + (A + fi 2 )R + (& + /ia)/? 2 , A = b - and /3 2 = b -@£&. 
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